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Abstract 

In the framework of the chiral quark-soHton model of the nucleon we investigate 
(NJ ', the properties of the polarized quark distribution. In particular we analyse the so 

' called anomalous difference between the representations for the quark distribution 

functions in terms of occupied and non-occupied quark states. By an explicit ana- 
<^ , lytical calculation it is shown that this anomaly is absent in the polarized isoscalar 

I distribution Au + Ad, which is ultraviolet finite. In the case of the polarized isovec- 

tor quark distribution Au — Ad the anomaly can be cancelled by a Pauli-Villars 



subtraction which is also needed for the regularization of the ultraviolet divergence. 



Oh! 

^ ; 1 Introduction 



Recently a rather successful program of computing the quark distribution functions in the 
framework of the effective quark-soliton model was developed ^, ^, |, |^ . The quark 
5^ i soliton model includes the chiral pion field U = and the quark field ip whose 

interaction is described by the Lagrangian 

L = i>{i-i^'d^- MW')i) . (1) 

In the mean field approximation (justified in the limit of the large number of quark colors 
Nc Pl) the nucleon arises as a soliton of the chiral field U 

U{x) = exp[i(n"r")P(r)], = r= \x\ . (2) 

This effective theory allows a quantum field-theoretical approach to the calculation of 
the quark and antiquark distributions in the nucleon. In contrast to naive quark composite 
models and to the bag model here we have a consistent approach reproducing the main 
features of the QCD parton model like positivity of the quark and antiquark distributions, 
various sum rules etc. 
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In terms of the quark degrees of freedom this picture of the nucleon corresponds to 
occupying with = 3 quarks the negative continuum levels as well as the valence level 
of the one-particle Dirac Hamiltonian H 

H = -i-i^-i^dk + M-i^W^ , (3) 

in the background soliton field U . For the pion field (Q) one can find the spectrum of the 
Hamiltonian 

H\n) = En\n) . (4) 

Various nucleon observables can be naturally represented as sums over eigenstates |n) 
of the Dirac Hamiltonian H. For example, the nucleon mass M^r is given by the sum over 
all occupied states or alternatively by the minus sum over all non-occupied states 

M^=^E„ = - Yl (5) 

n,occ n,non—occ 

In this expression the subtraction of similar sums is implied where the eigenstates \n) and 
the eigenvalues En are replaced by those of the free Hamiltonian, 

= |n^°^ ^0 = -^I'^l'^dk + M7°. (6) 

The physical reason for the existence of the two equivalent expressions in (|^) is that 
the polarized Dirac sea picture can be formulated either in terms of quark or in terms of 
antiquark states (occupied antiquark states correspond to non-occupied quark states). 

Formally the equivalence of two representations for Mjv follows from the identity 

^E„+ Yl E^ = YEn = TrH = 0. (7) 

n,occ n,non—occ n 

At the last step we took into account that the trace of H over the spin indices vanishes. 
Strictly speaking, this naive argument is not safe since the sums (H) over the occupied and 
non-occupied states are ultraviolet divergent and must be regularized. In principle, the 
ultraviolet regularization could lead to an anomalous difference between the summation 
over occupied and non-occupied states but in the case of the nucleon mass (||) one can 
check that in the regularizations like Pauli-Villars or proper-time ones the anomaly is 
absent. E.g. the proper-time regularized version of (|^ with the ultraviolet cutoff A is 

= 0. (8) 

We can reformulate this verbally as the "absence of the anomaly" in the nucleon mass 
M^v (in the proper time regularization). The usage of the word "anomaly" is invoked by 
the similarity with the axial anomaly which can be interpreted as nonvanishing trace of 
75 in e.g. the proper-time regularization. 

The main object of interest in this paper is the study of the quark distribution func- 
tions. In the mean field approach (justified in the large A^^^ limit) the quark distributions 
can be represented as single or double sums over occupied or non-occupied one-particle 
eigenstates (|^) of the Dirac Hamiltonian (^. We shall see that for the same parton distri- 
bution one can write two naively equivalent representations but whether this equivalence 



lim Y En exp -El/ A 



A — ».rv-^ 
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persists or not when one takes into account the ultraviolet regularization is a rather subtle 
question and the situation is different for different distributions. Moreover, even in the 
limit of the large cutoff, the cancellation of this anomalous difference between the naively 
equivalent representations is sensitive to the regularization used. 

Let us start from the unpolarized isosinglet quark distribution u{x) + d{x) which is 
given by the following expressions ([![]) in the leading order of the 1/Nc expansion 



n,nan—occ 



Also here the subtraction of similar sums with the eigenstates and eingenvalues of the 
Hamiltonian (Q) replaced by those of the free Hamiltonian is implied. The result (|^) 
has a transparent physical meaning of the probability to find a quark with momentum 
fraction x in the nucleon in the infinite momentum frame. In ref. it was shown that 
in the Pauli-Villars regularization the sums over occupied and non-occupied states in (|^) 
really give the same result. 

We stress that the fact of the equivalence of the two representations for parton dis- 
tributions is crucial for the positivity of unpolarized distributions and for the validity 
of various sum rules inherited by the model from QCD [|T|. Therefore the check of this 
equivalence is an essential part of the calculation of parton distributions in the chiral 
soliton model. 

Now let us turn to the polarized quark distributions. In the leading order of the 1/Nc 
expansion only the isovector polarized distribution survives 



Auix) - M{x) = -\n, y f 



(27r)3 \^ Mn I 

X (n|p)(l + 7V)^V(P|^) • (10) 

Compared to the expression (j^) for u[x) + d[x) here we have an extra factor which 
reflects the fact that now we deal with the isovector polarized distribution. The factor of 
1/3 comes from a careful treatment of the rotation of the soliton [|T|. 

One can ask whether the summation over the occupied quark states in (|T0|) can be 
replaced by the summation over non-occupied states 

1 Ar V- f d^p / p^ + En \ 

X ] 



Au{x) - Ad{x) = -N^ Y: Jj^.S 



n,non~occ 



(27r)3 \^ M, 



N 



X (n|p)(l + 7V)A^(p|^)- (11) 

In this paper we shall show that in the case of the Pauli-Villars regularization (the sum 
over states n in (^) is logarithmically divergent) the two representations (|TU|) and ( [TT| ) 
are really equivalent. 

We stress that the equivalence of the summation over the occupied and non-occupied 
states is very sensitive to the choice of the regularization. For example, if instead of the 
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Pauli-Villars regularization we simply cut the summation over quark states in ( ]T0| ) includ- 
ing only states with \En\ < ujq then a nonzero difference between the two representations 
(PUP and ( ]TT| ) will remain even in the limit of the infinite cutoff (X'o ^ oo. The mechanism 
how this anomalous difference appears is similar in many respects to the famous axial 
anomaly. In particular, such similarity manifests itself in the fact that the anomalous 



difference between the two representations (|TOD and ( pAf ) can be computed analytically in 



the limit ^ oo. The calculation of the anomalous difference is presented in this paper. 

Although the regularization including only states with \En\ < is not acceptable as 
a physical one and the Pauli-Villars regularization is more preferable in this respect, we 
want to emphasize that in the practical calculations based on the numerical diagonaliza- 
tion of the Dirac operator in the background soliton field, the \En\ < ujq regularization 
appears naturally. Indeed, in the numerical calculation one can work only with a finite 
amount of quark states so that one actually uses both Pauli-Villars subtraction (with the 
regulator mass Mpy) and the \En\ < Wq regularization. The pure Pauli-Villars subtraction 
is simulated by working with uo 3> Mpv- The numerical calculation is rather involved 
and the analytical result for the anomaly in the \En\ < ujq regularization is very helpful 
for the control of numerics even if the anomaly cancels after the Pauli-Villars subtraction. 

Now let us turn to the polarized isoscalar quark distribution Au{x) + Ad{x) which 
gets the first nonzero contribution only in the subleading order of the l/N^ expansion 

NrMN ^ ^ 1 



Auix) + Adix) = ^—^ y: E 



^-^ m,aU n,occ E^n 

X (ri|r^|m)(m|(l + 7V)7^^(^n + - xMjv)|n) 
+ § ^ E (^1(1 + lW)r'lH{E^ + P'- xMr,)\n) . (12) 

UX n,occ 

Here P"^ is the quark momentum projection on the third axis 

and / is the moment of inertia of the soliton 

Another representation for Au + Ad can be written in terms of the summation over 
non-occupied states n 

NMn ^ ^ 1 



Au{x) + Ad{x) = j: e 



m,aU n,non-occ Em 



X {n\T'\m){m\{l + -f^Y)rS{En + P- xMjv)|n) 
§^ E H(l + 7V)rV^(i5n + P'-xM^)|n). (14) 



n,non—occ 



The numerical calculation of Au + Ad with the Pauli-Villars subtraction was presented 
in paper [|]. Unfortunately there the question about the equivalence of the two repre- 
sentation (|12D and (|l|) was not investigated properly. Also the Pauli-Villars subtraction 
was used in paper H without proper justification. 
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In this paper we show that if one cuts the sum over occupied (non-occupied) states n 
allowing only < cuq in the eqs. (|12D, ([T^) then in the infinite cutoff limit uo ^ oo 

1) both representations (JT^) and (0) have a finite limit (i.e. A'u(x) + /S.d{x) has no 
ultraviolet divergences) , 

2) the two representations (|T2p, (|Tip give the same result. 

Comparing the last terms in the rhs of representations (0) and ([1^) for Am + Ad 
with expressions (|TUp and (|TT]) for An — Arf we see that the total expression for An + Ad 
contains a contribution proportional to ^ [An(x) — ^d{x)]. 

Therefore we start our analysis by investigating the anomaly of Am — Ad which we do 
in section |^. In section |^ we show by explicit calculation that for the quark distribution 
Am + Ad there is no anomalous difference between the summations over occupied and 
non-occupied states. In section ^ we discuss the numerical results and compare them the 
GRSV parametrization of experimental data. 



2 Anomaly of Aix(x) — A(i(x 

As it was explained in the introduction one of our aims is to investigate whether the two 
representations ( p!OD and ([l^) for the polarized isovector quark distribution Au{x) — Ad{x) 
are equivalent. The answer to this question is sensitive to the ultraviolet regularization. 
Let us start from the regularization that allows only the quark states n with \En\ < ujq- 
In this regularization eq. ([T0|) can be rewritten as follows. 



[Am(x) - Ad{x)] 



U)0 

\occ 



--N.Mr, 



dujTi 



5{H - uj)5{u + P3 _ a;Mjv)r3(l + 7V)75 



Here H is the Dirac Hamiltonian (^ and is momentum operator 
representation ([TT|) becomes 

[/\u{x) - M{x)Z[ 



lon—occ 



UJQ 



duTi 



5{H - uj)5{uj + P3 _ xMm)t^{1 + 7V)75 



(15) 
Similarly, 

(16) 



The main results of this section can be formulated as follows 



1) Both [Am(x) - M{x)]2. and [Au{x) - Ad{x 
in the hmit of large cutoff uq ^ oo 

[Auix) - Ad{x)]Z - [Am(x) - Adix)] 



LOO 

non—occ 



are logarithmically divergent 



N,MnM\ uo 

m — 

127r2 M 



d^k 
(2^ 



Sp.; (f/[k])+r^[/(k) e F - \x\Mj^) + 



(17) 
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2) In the difference [Au{x) — Ad{x)]'^°^ — [Au{x) — Ad{x)]'^l^_^^^ the ultraviolet di- 
vergences cancel and the Uq ^ oo limit of this difference reduces to the following finite 
expression 

l}m^{[Auix) - Adix)]Z - \^u{x) - Arf(x)]-„__} 



a;o- 
1 



127r^ 



log — V-V^^ — Spj, 



(27r 



AT 



*(^[k])+^(k) 



where ?7(k) is the Fourier transform of the chiral mean field ?7(r) entering the Dirac 
Hamiltonian fRT) 



f/(k) = y" d\e~'^^''^ \U{j) - 1] . 



(19) 



Note that \Au{x) — A(i(s)]^"^ and \Au{x) — Ad{x)Y^^^_^^^ separately are given by compli- 
cated functional traces (|15D and (|l^) which can be computed only numerically. The fact 
that the anomalous difference between the representations in terms of the occupied and 
non-occupied states reduces to a simple momentum integral (|1^) is highly nontrivial and 
is similar to the well known fact that the famous axial anomaly gets its contribution only 
from the simplest diagram. 

The fact that the divergence ( |17|) is proportional to means that this divergence 
can be removed by the Pauli-Villars subtraction so that the following combinations are 
finite 

\Au{x) - Ad{x)\Z 
- hm [Auix) - Ad{x)Cj' - [Au{x) - Arf(x)]-f 



lim <^ \Auix) - Adix)] 



[Au{x) - Ad{x)] 
M' 



uJo,M 
non—occ 



PV 

non—occ 
2 



Mly 



[Au{x) - Ad{x)] 



LL)o,Mpv 
non—occ 



(20) 



(21) 



Next since the anomaly (18) is proportional to we see that in the Pauli-Villars regu- 
larization the summation over occupied and non-occupied states gives the same results: 



[Au{x) - Ad{x)]Z = [^u{x) - Ad{x)] 



PV 

non—occ 



(22) 



Now let us turn to the derivation of the result (0) for the anomalous difference 
between the summation over occupied and non-occupied states. Subtracting (16) from 
(p!5|) we obtain 



i non—occ 



--N,Mn 



[Au{x) - Ad{x)]Z - [^u{x) - Ad{x)] 

6{H - uj)6{uj + p3 _ xM^yil + 7V)75 



dtoTi 



We use the following representation for the operator delta function 5{H 



UJ 



5{H 



signcj 
2Txi 



1 



(H + uj). 



(23) 



(24) 
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The squared Dirac Hamiltonian (^) is 

= -d^ + M'^ +iM{^''dkU^') . (25) 



Now (p3|) takes the form 



-U)0 



x5(cu + P3 - xM;v)r'(l + 7V)75} • (26) 

Next we expand the "propagator" in the rhs in powers of iM{'-j^dkU'^^) 

1 _ 1 

-(92 + M2 - - io + iM{-i^dkU^^) ~ -(92 + M2 - - io 

^ -.M(7^'9,f/-^)-— — + . . . . (27) 



-92 + M2 - tu2 - ZO " -(92 + M2 - Cj2 - iO 

The first nonvanishing contribution to (EO) comes from the term hnear in iM{'^^dkU"'^ 



9 ^° rl 

[An(x)-Arf(x)C-[A«(x)-Ad(a:)C_ = --iV,M;vIm / ^ 

— a;o 

x(c<; - i7°7''(9fc + 7°Mf/^5)(l + 7V)757'''5(^ + - xM;v)r^} • (28) 

Computing the trace over the spin indices and turning to the momentum representation 
according to (^9|) we arrive at 



7on—occ 



[Au{x) - M{x)]2, - [Au{x) - M{x)]Z 

= |iV,M^M2lm j ^fc3Sp^,{(f/[k])+[7(k)r3} 

r doj r d^p signer 6{u+p^ — xMn) 

^ J 2^ J (27r)3 |k + p|2 + M2 - u;2 _ ^0 |p|2 + - - iO ' ^ ' 

-UJO 

The calculation of the integral over uj and p is straightforward and in the limit t^o ^ 00 
one arrives at final result (|TB|). 

Note that the limit ujq ^ 00 should be taken after computing the u and p integrals. 
Otherwise if one first computes integrals u and at fixed p^ and with uq = 00 then one 
gets zero 

°° duj J dp^ signcj S{u! + p^ - xMn) _ , . 

2^ J 271 |k + p|2 + M2 - tu2 - iO |p|2 + M2 - cu2 - zO ~ ■ ^ ' 



-00 —00 
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Actually after integrations over oj (in the interval —uiq < oj < uiq) and the integration 
over Ip-*"! is restricted at large uoq to the interval 



2ujn min 



max IxM^vl, |xMjv + k 



(31) 



In the limit of large cutoff cuq — * oo this interval of p-*" is shifted to infinity which explains 
why the nonzero result ( |18| ) is compatible with the vanishing integral (|30|) . 

A similar phenomenon occurs with the p-*" integration region in the contributions 
coming from the higher terms of the expansion (^). However, since the integrands of 
these higher order terms decay faster at large Ip-*"! these higher order terms give vanishing 
contribution to the anomalous difference ([TSD in the limit of the large cutoff u;o ^ oo. 

Restricting the integration over uj in eq. ( [2B| ) to the interval —ujq < < or to 
< < ciJo we can investigate separate distribution functions [/S.u{x) — /S.d{x)\^l^ or 
[Am(x) — A(i(x)]^°„_Q^^. In this case one gets nonzero contributions from all terms of the 
infinite series (0). However, it is not difficult to check that only the first nonvanishing 
term of this expansion is logarithmically divergent in the limit of large cutoff uoq ^ oo and 
this logarithmic divergence is given by (0). This logarithmic divergence is proportional 
to and therefore in our previous calculation of Am — Ad we could regularize it by the 



Pauli-Villars subtraction. Moreover, since the anomaly 
is cancelled by the same Pauli-Villars subtraction pf. 



is also proportional to M it 



3 Cancellation of the anomaly of /\u[x) + Ad{x] 

Now we turn to the investigation of Au{x) + Ad{x). The cutoff version of (0) is 



N,Mn ^ ^ 1 



[Au{x) + Ad{x)]Z = ^ E E 



21 ^ ^ E — E 



X {n\r'\m){m\{l + 7V)7 (^(^n + P- xM^)|n) 
+ Tlic ^ (n|(l + 7V)rV5(K + P'-xM^)|r^). (32) 

-U)0<En<Elev 

Although we use notations corresponding to the discrete spectrum actually most of the 
spectrum is continuous. The singularities corresponding appearing to -C/m — E-^ are as- 
sumed to be regularized according to the principal value prescription. 
Making use of (|TB|) we find 

[Au{x) + Arf(x)]- = [Au{x) + Ad{x)t}r " ^ ^ [A«(x) - Ad{x)]Z , (33) 



where 

[Au{x) + Ad[x)]^r = ^Y. E 



X (n|r3|m)(m|(l + 7V)7^5(^n + P^-a;M7v)|n). (34) 
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Similarly ( |I^ ) leads to 

[Au{x) + Ad{x)] 



non—occ 



[Au{x) + Ad{x)] 



non—occ 



3 c) 



(35) 



where 



[Au{x) + Ad{x)] 



non—occ 



N.Mr, 



E E 



■"n 



X (n|r''|m)(m|(l + 7V)7 ^(^n + P - xMjv)|n) 



(36) 



We see that 



[Au{x) + Arf(x)]-», - [Au{x) + Arf(x)]2 
[A?i(x) + Ad{x)]^^^' - [Au{x) + Ad{x)] 



non—occ 



non—occ 



3 d 



AIMn dx 



{[Au{x) - Arf(x)]- - [A^(x) - Ad(x)]-„„_} . (37) 



Here 



[A..(a:) + Arf(x)]ir - + A^(^)]ii^-occ 



21 



occ 

E E 

rn —uio<En<uio 



-E'n ~ Em / py 



X {n\T-'\m){m\{l + 7V)7 '^(^n + P - xMjv)|n) . 



(3^ 



We remind that here the principal value prescription for {En — Em) ^ is implied. This 
can be rewritten in the form 



[Au{x) + Arf(x)]l'j"° - [Au{x) + Ad{x)] 



(l)'-'o 
non—occ 



xSp 



P.V. 



T^5{H -uj)+ 5{H - uj)r 



H -uo, 

X 5{uj + P3 _ xMn){1 + 7°7^)7^} 



4/ 



1 



duo 



H-uj 



P.V. 



4/ J 27r 



xSp 



1 



1 



if-a; + iO i^-w + zO 



(5(cj + - xMn){1 + 7V)7^ 



1^0 



+ 



iM^Nr f du 



4/ 



27r 



-a;o 



X Sp 



Hence 



H-Lj-tO H-u-iO 



[Auix) + Arf(x)]i'2"" - [Auix) + Adix)]l 



(39) 



(l)<^0 
non—occ 
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-Im- 



21 



X 



-Im- 



MmN^ 7 du 



J 27r yH'^ — uj^ — zOsignu; 

1 



21 



2'K 



signci^Sp 



X 



-Im 



i72 



M^Nr 7 dto 



21 



— UJ^ — lU J 



— u/u 

1 



-5(cu-z93-a;Mjv)(l + 7V)7' 



(40) 



-^2 + M2 - cj2 - iO + iM{-f''dkU^^- 

The rest of the calculation is similar to how we worked with expression (p6|) for the 
anomaly of Au{x) — Ad{x). 

Nonzero contributions to the anomaly come from the expansion of the propagators up 
to terms linear and quadratic in iM{'-)^d}JJ"'^): 



= Ai{x) +A2{x). 
Here Ai{x) corresponds to terms linear in iM{'~f^dkU'^^) 



l(l)i^O 



(41) 



Ai(x) = Im 



MnN^ 7 du 



271 



21 



signcuSp 



1 



-uio 



.Q2 + ^2 - a;2 - 



X (cj - i7°7'=9fc + -i^MW^)T^{uj - i7°7'=9fc + 7°Mf/^«) 
1 



X ■ 



-92 + M2 - c<;2 - 20 



iM{-f^diU"'') 



X 



-92 + M2 - CU2 - iO 

and is quadratic in iM^'j^dkU^^) 



5(cj-?a3 -xM^)(l + 7V)7' 



(42) 



T ^^^^ 7 • T / 1 



- Cj2 - iO 



2M(7™a^?7''^; 



1 



-92 + M2 - Cj2 - iO 
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-92 + M2 - Cj2 - zO 



1 



X (cj - i7°7'=9fc + 7°M[/^^)r=^(cu - i7°7'9, + 7°Mf/^'^0 
1 



X ■ 



-^M(7"9„x/^5; 



X- 



1 

-92 + M2 - CU2 - iO 
1 



iM(7'^a„f/^5 



X 



5(cu-za3-sM;v)(l+7°7')7' 



-92 + M2 - cj2 - iO 
A straightforward calculation leads to the following results for Ai{x) and A2{x) 



Aiix) 



A2{x) 



(27r)3P 



In 



1 + 



MnNM'^ f d^k 



87r2/ 



i2n) 



xMn 
In 



Sp [f/(k) ' f/(k) 

+ xMn 



xM 



N 



I Id 

X T^ + -- 



Now we insert these results into (O) 



Sp [/(k) f7(k) . 



[Au{x) + Aci(x)]^r - + M{x)t^_^^ 



MnN.M^ r d^k 
167r2/ J (2^ 



In 



(l)<-'0 
occ 

+ xMtv 



l(l)(^0 



xM. 



N 



_d_ 

dk^' 



+ ~ 



Sp f/(k) f/(k) 



Note that shifting the integration variable 

k^ k^ - xMn 

we obtain 



In 



d^k 
(2^ 



1 d r d^k 



k^ + xM, 



N 



_d_ 



Sp Uik) U{k) 



MnOxJ {2tt 



In 



P + xMj^ Sp f7(k) t/(k) . 



Therefore 



[Auix) + Arf(x)]ii^ - [Au{x) + Arf(x)]^L„__ 



Ar,M2 d r d^k 



IQtt^I dx J (27r 



In 



k^ + xM, 



N 



xMn 
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(l)'-'o 
non—c 

Spjr^ [t/(k)l^f7(k)l . 



Inserting this result and (|I8|) into ( p7| ) we observe a complete cancellation: 

[Au{x) + Ad{x)]Z - [Au{x) + Arf(x)]-„__ = (50) 

Thus the isoscalar polarized quark distribution Au{x) + Ad{x) is nonanomalous. 

Using similar methods one can check that function Au{x) +Ad{x) is free of ultraviolet 
divergences: although the two separate terms in the rhs of (^3]) are UV divergent the total 
sum is finite. 



4 Numerical results 

The numerical results for the isovector polarized distribution function Au{x) — Ad{x) 
are given in [Q. For the computation of Au{x) + Ad{x) (|T^, (|T^ we use the numerical 
methods which were developed in and later extended in for the computation of the 
isovector unpolarized distribution. 

The eigenvectors and eigenvalues of the Dirac Hamiltonian (|^) are determined by 
diagonalizing in the free Hamiltonian basis (|]) . This basis is made discrete by placing the 
soliton in a three-dimensional spherical box of finite radius D and imposing the Kahana- 
Ripka boundary conditions Both Au{x) — Ad{x) and Au{x) + Ad{x) were computed 
using the standard value of the constituent quark mass M = 350 MeV as derived from 



the instanton vacuum 12 



In our calculation we use the self-consistent solitonic profile P{r) (see e.g. ref. IH 



However, performing the numerical calculations in the finite spherical box one should be 
careful about the large distance effects. To be safe, we artificially exponentially suppress 
the pion tail of the soliton profile at large distances so that the field vanishes outside the 
box (a similar problem in the calculation of qa was studied in |ll5| ). 



In Fig. |1| we compare our numerical results for the anomaly of Au{x) — Ad{x) with 
the analytical result (p^ . We observe a rather good agreement. 

Fig. shows the numerical results for the Dirac sea contribution to Au{x) + Ad{x) 
based on the two representations (occupied and non-occupied). We see a reasonable agree- 
ment between the two results which confirms the absence of the anomaly in Au{x)+Ad{x). 
Some difference between the two curves at negative x is finite-box artefact. Increasing 
the size of the box one can see that this difference tends to disappear. 

In Fig. 3 we compare the result of the calculation of Au{x) + Ad{x), Au{x) + Ad{x) 
with the GRSV-LO parametrization |jl8[| at the low scale of the model /i = 600 MeV. We 
see that the quark distribution Au{x)+Ad{x) is in a reasonable agreement with the GRSV 
parametrization whereas the antiquark distribution Au[x) + Ad{x) obtained in the model 
is considerably smaller than that of the GRSV parametrization. Note that the polarized 
antiquark distributions are not directly accessible in inclusive hard reactions. Due to 
the lack of data the GRSV parametrizations therefore are based on certain assumptions, 
e.g. in the GRSV analysis it was assumed that Au{x) = Ad{x). In contrast to this the 
QCD large Nc counting and the quark soliton model predict a large flavour asymmetry 
Am(x) > Ad{x) in the light polarized sea. Some physical applications of this have been 
studied in refs. 0, HI, H. 
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Fig. 4 shows our predictions for the polarized antiquark distributions Au{x) and 
Ad{x) separately at the scale fi = 600 MeV. 

Since the quark distribution Am + Ad is finite, no ultraviolet regularization is needed 
for this quantity. There is even an argument against regularizing Au + Ad coming from 
the fact that the first moment of this distribution is related to the imaginary part of the 
quark determinant in the background soliton field which has to be left nonregularized if 
one wants to keep baryon number conserved - this is an analog of the nonrenormalizability 
of the Wess-Zumino term in pure chiral models. 

Several comments should be made about the calculations of Am + Ad within the 
same model by Wakamatsu et al. who published three different versions of the calcula- 
tion in papers ^, [1^. In paper [|l^ one of the terms was overlooked. This mistake 
was corrected by the authors of 0]. The revised version of calculation of Wakamatsu 
et al. was published in p. In this paper the question about the anomalous difference 
[Au{x) + A(i(x)]^°^ — [Au{x) + Ad{x)]'^l^_^^^ was investigated only numerically but the 
accuracy of the calculation did not allow the authors to draw any conclusions concerning 
whether this difference vanishes or not. Actually the numerical accuracy of the agree- 
ment between the two representations which we observe in our calculation (see Fig. Q), 
and which is necessary for a proper evaluation of the parton distributions, is of two or- 
ders magnitude better than the same difference presented in 0. The practical solution 
accepted in ^ was to use [Au{x) + Ad{x)]'^°^ for x > and [Au{x) + Ad{x)]'^l^_^^^ for 
X < (i.e. for the antiquark distribution). As it was explained above, Au{x) + Ad{x) 
should not be regularized contrary to what the authors of ||^ do. 

The polarized structure functions were also estimated in the work [^] in the Nambu- 
Jona-Lasinio model within the valence approximation. In ref. [|l| it was shown that the 
valence approximation leads to a number of inconsistencies: antiquark distributions are 
negative, sum rules are violated, etc. 

Our work has been performed within the quark-soliton model with two quark flavors. 
In the case of the model in the flavor SU (3) the same quantity should be interpreted as 
Au + Ad + As. 

The flrst moment of the Au{x) + Ad{x) gives the singlet axial charge. Our result of 
9a^ = j\ dx{Au + Ad){x) = 0.35 agrees with the calculation performed in other works 



10, Note that in the calculation of this charge no ultraviolet regularization was used. 



5 Conclusions 

We have proved that the representation of singlet polarized (anti) quark distributions in 
the chiral quark-soliton model as a sum over quark orbitals is ultraviolet flnite and free of 
quantum anomalies. This is a serious check of the consistency of the quark-soliton model. 

In fact, the cancellation of quantum anomalies in the model is related to the fact 
that certain basic properties of QCD as a local quantum fleld theory are realized in 
the model. The equivalence of the summation over occupied and non-occupied states is 
directly connected to anticommutativity of fermion flelds at space-like intervals. Actually 
this locality property has direct relation to the positivity of quark and antiquark densities 
in the quark soliton model [jl], ||]. 
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Another consequence of the cancellation of anomalies is that the model results for the 
parton distributions are compatible with the charge conjugation invariance: the quark 
distributions in nucleon coincides with the antiquark distributions in the antinucleon. 

From the practical point of view the results presented in this paper allow us to conclude 
that for the calculation of the singlet polarized quark and antiquark distributions no Pauli- 
Villars subtraction is needed. Additionally the numerical check of the cancellation of the 
anomalies is a powerful tool to control the accuracy of the numerics. 

We have computed the singlet polarized quark and antiquark distributions which arise 
in the subleading order of 1/iVc expansion. We found the quark distribution Au{x)+Ad{x) 
to be in a reasonable agreement with GRSV |]T8|] parametrization of parton distributions 
at low normalization point. A remarkable prediction of our model is that the polar- 
ized distributions of u and d antiquarks are essentially different, see Fig.4. Usually, in 
parametrizations of polarized parton distributions, it was assumed that Au{x) = Ad{x), 
which is not confirmed by our model calculations (see Fig. 4). It would be extremely 
interesting to include into the fits of the data the flavour decomposition pattern for po- 
larized antiquarks obtained in our model calculations. Future experiments at HERA and 
RHIC investigating Drell-Yan lepton pair production in polarized nucleon-nucleon colli- 
sions will clarify the situation. For a discussion see |jT9|, |20| . Let us note that in the singlet 
polarized channel under the evolution the quark distributions mix with polarized gluon 
distribution. Analysis of ref. [^, |l| in the framework of the instanton model of the 
QCD vacuum shows that the gluon distribution is parametrically smaller (suppressed by 
M'^/Mpy) than quark and antiquark distributions. In order to obtain a non-zero result 
one has to go beyond the zero-mode approximation of ref. |]23| and/or consider contribu- 
tions of many instantons. Both ways would lead to extra powers of the packing fraction of 
instantons. This means that gluons at low normalization point inside the nucleons appear 
only at the level of M'^/Mj.y. 
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Figure 1: Analytical (solid) and numerical (dashed) results for the anomalous difference 
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Figure 2: Results for continuum contribution [Au + Ad]sea based on the occupied and 
non-occupied representations. 
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Figure 3: The quark soliton model results for a;[AM-|- Ac?] and a;[Ati-|- Ad] versus LO-GRSV 
parametrization at the scale jj, ~ GOOMeV. 
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Figure 4: The quark soliton model predictions for xAu and xAd at the scale /i ~ GOOMeV. 
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